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Abstract 

The nuclear reaction e + 3 He + 4 He — > e + 7 Be is considered at thermonuclear 
energies. The motion of the electron is treated within the adiabatic approx- 
imation and the 3 He — 4 He scattering state is calculated using an effective 
3 He — 4 He potential constructed via the Marchenko inverse scattering method. 
The reaction rate thus obtained for solar interior conditions is approximately 
10 -4 of the corresponding rate for the radiative capture 3 He + 4 He — > 7 Be + 7. 
However, at the conditions corresponding to the early stages of the primordial 
nucleosynthesis, the non-radiative fusion is five times faster than the radiative 
one. The importance of the non-radiative processes for the energy balance in 
the sun, which for electrons is totally opaque, is discussed. 

PACS number: 21.45.+V, 95.30.-k, 96.60.-j, 97.10.Cv 

Key words: Thermonuclear reactions, fusion, stellar plasma, primordial nu- 
cleosynthesis, pp-cycle. 



I. INTRODUCTION 



In the standard model of the sun it is assumed that all nuclear reactions forming the pp- 
chain, are generated by various two-body nucleus-nucleus and electron-nucleus collisions 
with only one exception, namely, the e + p + p — > d + v reaction [JTJ] . Other three-body 
reactions are ignored on the grounds that the triple collisions occur less frequently than the 
binary collisions. 

Defining the collision as a situation in which the particles approach each other to a 
distance less than the range R of their interaction, then the frequency of two-body collisions 
per unit volume is 

T\i = ■KR 2 vi 2 n 1 n 2 , 

where V\2 is the relative velocity of particles 1 and 2, and Uj are their densities. The frequency 
of the triple collision JF 12 3 is obtained if we multiply JF 12 by the probability of finding the 
particle 3 present when 1 and 2 are scatterred ||, viz. by (Atc /3)R 3 n 3 . Then 
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4tt 2 , 

From these elementary considerations one can find that triple collisions are less probable 
than the binary ones. For example, in the centre of the sun the densities of electrons, 3 He, 
and 4 He nuclei are (these can be derived from the Tables 4.4 and 4.5 of Ref. [T|]) 

n e = 0.60 x 10 26 cm" 3 , 
n 3Hc = 0.23 x 10 21 cm" 3 , (1) 
ri4 He = 0.14 x 10 26 cm" 3 , 

which together with a nuclear reaction range of R ~ 10 fm and a velocity of t>3He 4 He ~ 
3 x 10 7 cm/sec corresponding to the temperature 1.56 x 10 7 ° K, give 

-^3 He 4 He ~ 0.3 x 10 30 cm- 3 sec _1 , 

He 4 He ~ 0.8 X 10 2 ° CHI 3 SeC 1 . 

This, however, is not a sufficient proof that the three-body reactions can be ignored. This is 
due to the fact that not every collision results in a nuclear transition. The actual reaction rate 
depends also on the transition probability which may be higher for three-body rather than 
for two-body initial states. This stems from the fact that triple collisions are kinematically 
less restricted than the binary ones. Some binary reactions are forbidden by conservation 
laws, e.g., conservation of angular momentum, parity, isospin etc. However, such reactions 
could take place in the presence of a third particle. 

Calculations of the three-body reaction rates performed so far |3|-|5| , corroborate this by 
showing that the simplified considerations based on the frequency of collisions overestimate 
the differences between the two- and three-body reaction rates. Thus, for example, in Ref. 
|| it was shown that under the conditions which prevailed during the primordial nucleosyn- 
thesis, the non-radiative capture e + n + p — > e + d occurs at a reaction rate which is five 
orders of magnitude less than the corresponding rate for the radiative capture n+p — > 7+d. 
In Refs. j|,[5| the reaction rates for the non-radiative processes e+p + d—>e + 3 He and 
e + p + 7 Be — ► e + 8 B were calculated. Again, it was found that they are only four orders of 
magnitude less than for corresponding binary reactions p+d — > 7+ 3 He and p+ 7 Be —>■ 7+ 8 B. 
In contrast, a simplified estimation provides a rate which differs by ten orders of magnitude 
from them. 

The above examples clearly demonstrate that a reliable answer to the question of whether 
or not a particular three-body reaction is negligible, can be obtained only when proper 
calculations are made. And surely in the case under consideration, such calculations should 
be done as we have the 7 Be paradox on top of the long-standing solar neutrino problem. 
This paradox originated from the combined analysis of all experiments, in which the neutrino 
flux from the sun was measured. This analysis came to the conclusion that the production 
of 7 Be nuclei (more precisely, the flux of neutrinos due to 7 Be reactions) must be strongly 
suppressed or even be negative ffJJ^]. This implies that something is wrong either in the 
standard model or in the experimental data. In this respect the fate of 7 Be in the pp-chain 
is of special interest. 

According to the standard model this nucleus is produced via the radiative fusion 
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3 He + 4 He — > 7 Be + 7 , (2) 

in which the released energy (1.59 Me V) is carried away by the photon. Another possibility 
for a 7 Be production is the three-body non-radiative fusion 

e + 3 He + 4 He — ► e + 7 Be . (3) 

This is a kind of Auger transition in which the same amount of the released energy is carried 
away by the electron. 

As we have seen, the triple collisions (^) are less frequent than the binary ones (^) and 
thus one can argue that they are unimportant. We have, however, at least three reasons 
to perform explicit calculations for the reaction rate of the process Firstly, since the 
neutrino paradox persists we have to be scrupulous in examining any conjecture on the 
matter. Secondly, although the three-body reaction may be unimportant in the neutrino 
problem, it may contribute to the energy balance in the sun since the energy released in 
each non-radiative fusion event remains in the center of the sun because the electron which 
carries it, cannot escape. In contrast to photons, for electrons the sun is totally opaque. To 
estimate this additional heating, we need to know the three-body reaction rate. Finally, at 
the early stages of the primordial nucleosynthesis the temperature as well as the particle 
densities were much higher than they are in the solar plasma. We therefore may expect that 
the three-body reactions were important at that times. 

The paper is organized as follows. In Sec. II we describe our formalism in general and in 
Sec. Ill we outline the procedure and approximations employed to evaluate the transition 
operator. Sec. IV is devoted to the potentials we use to describe the interactions among the 
particles, and In Sec. V we present our results and conclusions. 

II. THREE BODY REACTION RATE 

The nuclei 3 He and 4 He are stable clusters (in what follows we denote them by r and 
a respectively) in the low-lying states of 7 Be (see Ref. || and references therein). We 
can, therefore, reduce the eight-body problem e(NNN)(NNNN) to the three-body one, 
namely, e — r — a. The corresponding Jacobi vectors in configuration and momentum spaces 
are shown in Fig. [TJ 

Let 7£(k, p — > p') be the reaction rate per unit volume per second for the collision process 
(0) which starts with the momenta (k, p) and ends up with p'. The plasma particles are in 
a thermodynamical equilibrium and their momenta are distributed according to Maxwell's 
law 

iV k (0) = (2nfiK6)- 3/2 exp 

N p (9) = (2imiK6)- 3/2 exp 

where and N p are the probability densities, \i is the ra-reduced mass, m is the electron 
mass, k is the Boltzmann constant, and 9 is the plasma temperature. We are cocerned with 
the total rate of the transition from an initial state with any (k, p) to a final state with all 
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possible p'. Thus the reaction rate 7Z(k, p — > p') must be averaged over the initial momenta 
k and p and integrated over the final momentum p', i.e., 

(n) = J J jdkdpdp'n(k,p^p')N k (e)N p {9). (4) 

Similarly to the two-body reaction theory where the average reaction rate (lZi2)e is written 
as a product of (012^12)0 (which is referred to as the reaction rate per particle pair) and the 
particle densities ||, 

(Ki2)e = "1^2(0-12^12)0 , 
the three-body reaction rate (TZ)e can also be factorized in the same manner, 

(TZ) e = n e n T n a (J:) e . (5) 

Indeed, using the general formula for a three-body reaction rate M we obtain 



K (k, p - p') = ^n e n T n a V 5 - E^ - — - ^ 
1 J 2 \^ mj V 2 ™ 2/i 2m, 

X I i^Tjm^P' \ T \ ^ m s]P)\ 2 , 



(6) 



where the 5-function secures the energy conservation for the transition to the n-th level of 
7 Be nucleus (see Fig. 0) accompanied by the energy release and T is the transition 
operator for the process (Q). The initial and final asymptotic state wave functions for 
this reaction, |k, m s ;p) and |V4jm) P'K are characterized, apart from the momenta, by the 
nuclear spin s = 1/2 of the 3 He, the total angular momentum j of 7 Be (which depends on 
the state n), and their third components m s and rrij respectively. They are normalized as 

( k', m' s ; p' | k, m s ; p) = 5 m>ls 5(k - k') 8(p - p') , 



( V'SioP' I ^7iL,;P) = 8 n'n 8 m >. m . 8{p- p) . 

The sum | Z)m s m- m ® * s f° r averaging over the initial and summing over the final spin 
orientations. Therefore the quantity (£)# defined by Eq. @, in analogy to the two-body 
case, may be referred to as the reaction rate per particle trio and is given by 

( S ) e = i^l! Y, I I [ dkdpdp' 8 - E^ n 

2 nm 3 mj J ^ •* \2lTl 

x I (4%.; p' | T | k, m s ; p) | 2 N k (6)N p (9) . 

In the next two sections we describe how the various ingredients needed to obtain the matrix 
element (ip7j mj ', p' | T | k, m s ; p) were calculated. 
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(7) 
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III. TRANSITION OPERATOR 



The three-body quantum state |e + 3 He + 4 He) belongs to the Hamiltonian 

H = H + h + V N + V e , (8) 

where H and ho are the kinetic energy operators associated with the Jacobi variables r and 
p respectively (see Fig. [[J), 

v N = v {s) + v {c) 

JM "rot ' to 

is the nuclear t-ol potential, which includes strong and Coulombic parts. The V e consists of 
the e-r and e-a Coulomb potentials 

V e = V eT + V ea . 

The various transitions in this three-body system are determined by the relevant matrix 
elements of the T-operator obeying the Lippmann-Schwinger equation 

T{z) = V N + V e + (V N + V e ) - 1 -T(z) . (9) 

z — tt — n 



The probability of the transition (|3|), i.e. of 



|k,m s ;p) W ^;p) , (10) 

is therefore determined from the matrix element (V4jm-j P' I ^ I k, m s ; p). 

Making use of the special conditions prevailing in solar plasma and the smallness of the 
electron mass, we can significantly simplify Eq. (|J). Indeed, the average kinetic energy of 
the particles in the plasma, < E km >~ k6 w IkeV, is the same for nuclei and electrons, 
but the velocity of an electron is three orders of magnitude higher than that of 3 He or 4 He. 
Therefore, while the r approaches a very slowly, the electron dashes nearby picking up the 
energy and leaving the heavy particles in a bound state. Due to this we can treat the relative 
ra motion adiabatically. For a three-body collision such an approximation results in two 
simplifications. In the first one the action of the potentials Vn and V e can be separated. 
Indeed, while the electron starts from its asymptotic state |p), the heavy particles have 
already interacted via the potential Vn and formed the two-body scattering state 

|k,m s ) ==> |-0k,m s ) • (11) 

Therefore, instead of the transition ( |10D caused by both Vn and V e , in the adiabatic approx- 
imation we may consider the transition 

IVv, s ;p> ^ 141,; p') , (12) 



where the interaction Vn is taken into account by flTTp. In other words, the transition ( |T0| ) 
effectively occurs in two steps, (JTTj) and ([12]). As a result the two-body scattering problem 
(|TT|) can be solved separately. 

In the second simplification the transition (|T2| ) may be described using the fixed scatterer 
T-matrix defined as 
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f(z) = V e + V e —^—f(z) . (13) 
z — n 

Within the above approximation we have 

(^7jmj;p'\T\k,m s ; p) « p'\f \ip Kms ; p) . (14) 

A further simplification can be achieved when Eq. ( |13"D is solved iteratively. For the solar 
plasma electrons one has the condition e 2 /hv < 1, which is a sufficient to treat Coulomb 
interactions in Eq. (O) perturbatively, i.e., 



f(z) = V e + V e —^—V e + V e —^—V e —^—V e + ■ • ■ , 

Z — tlQ Z — lift Z — iIq 

Furthermore, the average potential energy of the electron is of atomic order of magnitude, 

< V e >~ 10 eV, while its kinetic energy in solar plasma is two orders of magnitude higher, 

< E^ m >~ 1 keV, which implies that the above iterations should converge very fast. There- 



fore, we may retain only the first (Born) term (10 



Finally, the exact matrix elements of the T-operator can be replaced by the approximate 



ones 

',(") . JlTll, ™ . ~\ ~ /„/.( n ) 



{Wjmii p' I T I k, m s ; p) w mjmpP' I V e I ip^ ms \ p) • (15) 

This can be considered as a three-body generalization of the distorted wave Born approxi- 
mation (DWBA) which is widely used in the theory of the two-body nuclear reactions. 



IV. POTENTIALS 

The evaluation of the matrix element fll5|) requires the knowledge of the potential V e and 
the wave functions ipi^ ma (r) and tpTj m . (r) . In order to obtain the wave functions we solved 



the toi two-body problem using the Jost function method, proposed in Refs. [p^lg] , with 
the rot effective potential V TCe = V$$ + V$. 



A. Coulomb forces 

The electron as a point-like particle generates a pure Coulomb field. However, the electric 
charges of r and a are distributed within their nuclear volumes which have a size comparable 
to the typical range of the nuclear reaction (few fm). Furthermore, the average collision 
energy, ~ 1 keV, is comparable to the height of the Coulomb barriers and thus the charge 
distributions should be taken into account. This is achieved by assuming that the nuclear 
charge distributions are spherically symmetric and have a Gaussian form 

Qi(Q) = 2e (^y) * exp (-^) , i = r, a , (16) 
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where g is the distance from the centre of the nucleus. The normalization coefficient and 
the slope parameter in ([H]) are chosen to give the total charges and the r.m.s. radii of the 
nuclei, i. e. 



1 

2^ 



J qi(g) d 3 g = 2e , 
g 2 qi (g) d 3 g=(R 2 } 



For the r.m.s. radii we used the experimental values |TT|JT^1 

1.67 fm . 



'<J$) = 1.93 fin, y/(E*) 
Using the space charge distribution fll6|) , we obtain 



2e 2 / 
V ei {g) = — erf I ^ 



r, a 



(17) 



v%Kq) 
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n(Rl) Jo 



exp 



3(g' 2 + Q 2 ) 
2<J^> 



sinh 



gg 



erf 




. (18) 



At large distances this electron-nucleus potential (|T7j ) has the usual Coulomb tail 2e 2 /g. At 
short distances, however, it is essentially different without having a singularity at g = 0. 
The same can be said about the nucleus-nucleus interaction flTg|). 

The electron which participates in the three-body reaction, is only one of the great 
number of the plasma electrons. All the others are spectators surrounding the nuclei and thus 
their combined Coulomb field reduces the electric fields of the nuclei. In the standard Debye- 
Hiickel theory |l3j such screening effect is taken into account by introducing an exponentially 
decaying factor into the Coulomb part of the nucleus-nucleus potential. Following this 
approach, we replace the potentials (|T7|) and (|18|) by the screened ones 



V ei (g) V ei ( g) exp ( - 



exp \ - — 



r, a 



(19) 
(20) 



with the Debye radius being D = 21800 fm which corresponds to the solar plasma conditions 
and is typical for other stars as well JT5 . 



B. Nuclear forces 

As can be seen in Fig. ||, when 3 He and 4 He nuclei are fused in solar plasma, they may 
form either the ground or the first excited state of 7 Be, the quantum numbers j n being 3/2~ 
and l/2~ respectively. The rest of the excited states are situated very high and hence we 
can safely ignore the virtual transitions via them. Therefore, in constructing the nucleus- 
nucleus scattering wave function ^k,m s we assume that 3 He and 4 He may interact if their 
total angular momentum j < 3/2. This restricts the allowable values of their relative orbital 
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momentum t to t < 2. The state with £ = 2 can also be ignored since at collision energies 
~ IkeV the contribution from the higher partial waves diminishes very fast. Thus for the 
nuclear forces between 3 He and 4 He we consider those corresponding to the lowest three 
(£, j)-states, namely, (0,1/2), (1,1/2), and (1,3/2). 

To construct the corresponding ra-potentials, we employ the Marchenko inverse scatter- 
ing method []T7| , p~8|] which is briefly outlined below. Within this method we obtained energy 
independent local potentials V^(r) for each of the above three partial waves, which repro- 



duce the the Resonating Group Model (RGM) ra-scattering phase-shifts [|19|,^(J anc ^ §i ve 
the correct binding energies for the ground and first excited states of 7 Be (in the ra-model). 

In the Marchenko inverse scattering method a unique, energy-independent, ^-dependent, 
local potential Ve(r) is obtained from 

V e (r) = -2^-K e (r,r) (21) 
dr 

where the kernel Ki(r,r') obeys the Marchenko fundamental equation 

K e (r, r') + F e (r, r') + / K e (r, r")F e (r", r')dr" = 0. (22) 

Jr 

The driving term Fi(r,r') is given by 

F e (r,r') = — / wj(kr)[l - S e (k)}w£(kr')dk + A e wf (b t r)wt (b e r') . (23) 



2tt 



oo 



St(k) is the S-matrix for the specific partial wave, and the function wf(z) is related to the 
spherical Hankel function h\ (z) by 

wj(z)=ie M zh { +\z). (24) 



Furthermore, is the so-called asymptotic bound state normalisation constant, while b 



■2nE^, El being the ra bound state energy and \x their reduced mass. 
The evaluation of Fi(r, r') requires the knowledge of the S-matrix for all real energies from 
the elastic scattering threshold to infinity, together with the reflection property S^—k) = 
1/St(k), as well as the binding energies and the corresponding asymptotic bound state 
normalisation constants. It is greatly simplified by choosing a rational (Bargmann-type) 
parametrisation 

n ; ^k-ibf^k-ai 1 J 

where A^^ is the number of bound states present. The number Ni + Nb must be even 
to satisfy the properties of the 5-matrix. The a l n are complex numbers used to fit the 
(numerically) given S-matrix. We mention here that the so-called Pauli Forbidden States 
(PFS) are ignored. This implies that the resulting potentials are the unique supersymmetric 
shallow [[Zl]] partners sustaining only the physical bound states. 



With the choice (|25|) the integration in Eq. (p3|) can be easily performed analytically 



(for more details see |T8| , |2T 



The potentials thus obtained for the ra system in the channels Si/ 2 , P1/2, and P 3 / 2 
are shown in Fig. E| 
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V. RESULTS AND CONCLUSIONS 



The plasma temperature in the center of the sun is 9q = 15.6 x 10 6 °K [jTJ and decreases, 
by a factor of ten towards the surface. Calculating the reaction rate for the process @ 
we need, therefore, to cover this range of 9. The same pp-chain reactions occur in other 
stars as well with similar conditions. In the blue stars, which are hotter than the sun, 
the plasma temperature is much higher. Thus, it is worthwhile to extend our calculations 
beyond the solar temperatures too. Another reason for such an extension comes from the 
necessity to examine the role of the triple collisions in the primordial nucleosynthesis which 
also proceeded via the pp-cycle. The temperature at which the creation of the nuclei occured 
in the early universe was ~ 10 9o K ||. 

We, therefore, calculated the three-body reaction rate for the values of 9 ranging from 
1 x 10 6 °K to 3 x 10 9 °K. The results of our calculations are given in Table |I[ Since the 
reaction rate obtained can be used not only in models of the sun but also for other stellar 
objects as well where the plasma densities may be quite different, we present it in units of 
cm 6 mole _2 sec _1 which are convenient for a general use |23]. The meaning of these units is 



that instead of n e n T n a we multiply (£) by N\, i.e. the Avogadro number squared. Then for 
any specific densities the reaction rate, in the units cm _3 sec~ 1 , can be obtained from Table 
| by multiplying the values given there, by n e n T n a /N\. 

To answer the question of the role played by the three-body reaction (|3|) in the solar 
pp-cycle, we need to compare its reaction rate with the corresponding rate of the two-body 
process (0). Using the particle densities ([!]), we obtain for the center of the sun 

(ReTa}e = 7.15 x 10 3 cm" 3 sec _1 , 

while for the corresponding two-body process 



Their ratio is rather small, 



(^ra)flo = 2M x 107 cm- 3 sec _1 . 



(R-TOt)t 



2.94 x 10~ 4 , (26) 



but not as small as one would guess from the simplified considerations discussed in the 
introduction. The reaction rate {lZ era )g in other parts of the sun together with its ratio 
to the two-body reaction rate are given in Table [TT]. In obtaining these data, we used the 
dependences of n e , n T , n a , and 9 on the distace R from the solar centre given in Ref. Q. 

When we move away from the center of the sun the temperature and the electron density 
decrease. As a result the ratio (TZ eTa )g/(TZ Ta )g becomes even smaller than (p6|). But the 
nuclear burning occurs mainly in the central zone of the sun (R ~ 0.1R & ). For example, the 
two-body reaction rate (TZ TOl )g at the radius 0.1R Q is already one order of magnitude less 
than at the centre ||23|| . Meanwhile, within the distance R < 0.1R Q from the solar centre 
the ratio (TZ eTa ) g / ' {TZ Ta ) g remains practically the same, ~ 10~ 4 . Therefore, in calculating 
the network of the solar pp-chain reactions one can ignore the non-radiative synthesis of 7 Be 
only if the required accuracy is less than ~ 0.01 %. 

For more hot and dense stars such accuracy restrictions are more essential as the ratio 
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\ = 7 r We ( 2T ) 

is proportional to the electron density n e . Furthermore, by comparing the data of Table | 



with the corresponding data for the two-body reaction given in Ref . [£3|] , one can find that 
the ratio (S eTQ ,) e / (cr Ta v Ta )g increases with temperature. 

Another theory in which the omission of the three-body non-radiative synthesis of Be 
may be detrimental is that of the primordial nucleosynthesis during which practically all the 
nuclei were created when the temperature was in the range 

0.3 x 10 9 °K < 9 < 3 x 10 9 °K (28) 

with the electron density being (see Ref. f3j) 

9Q / 9 \* ( 5.93 x 10 9 \ , , . 

n e ^1.5xl0 29 [— J expl J cm" 3 . (29) 

At 9 = 10 9 °K, for example, formula (25) gives 



n e « 0.40 x 10 27 cm" 3 , 9 = 10 9 °K , 

which is one order of magnitude higher than the electron density ([TJ) in the solar interior. 

Using the data of Table | and the corresponding data for the two-body reaction (0) 
of Ref. |23[| , together with the electron density (p9|) , we calculated the ratio (p7|) for the 
primordial nucleosynthesis. The temperature dependence of this ratio which we found for 
the temperature range (p8|) is given in Table [ITJ. ^ * s seen that at the earliest stages of the 
nucleosynthesis, when the temperature was 9 « 3 x 10 9 °K, the non-radiative synthesis (|3|) of 
7 Be nuclei was five times faster than for the two-body fusion (0). Thus, if the electron density 
given by formula ( p9|) is not far from the actual density that existed at the early stages of the 
Universe, then three-body reactions of the type fl3|) must be taken into account in theories 
concerning primordial nucleosynthesis. It can make a differences since the uncertainty in 
the 7 Be production in the primordial nucleosynthesis is currently assumed to be of the order 

ofi6%m. 

Despite the fact that in the solar plasma the contribution of the three-body process is 
rather small, in enegy balance considerations this might not be the case. In each event (|3|) 
or (0) an energy of 1.59 MeV is released. The energy produced by the radiative capture (0), 
is carried away by a photon which can penetrate the inner layers of the sun and even escape. 
In contrast, the electron which receives the energy produced in the reaction (|3|) remains in 
the interior of the sun and eventually distributes its excess energy among the other plasma 
particles. This causes an additional heating which is not taken into account by the standard 
model of the sun. 

To evaluate the importance of three-body reactions in the energy balance quantitatively 
we calculated the rates of the energy released, due to the reactions ([3]) per cm 3 per second 
at different radial positions in the sun. These are presented in Table ||. These results were 
obtained using the reaction rates of Table [I] together with the radial dependences of the 
particle densities and the temperature given in Tables 4.4 and 4.5 of Ref. |TJ. By integrating 
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over the sun volume, we obtain the total energy generated in the sun every second from the 
non-radiative synthesis of 7 Be, 

( = 0.582 x 10 16 ^ 

V at I sec 

V / era 

As compared to the total photon luminosity of the sun Q 

L & = 0.386 x 10 34 — 
sec 

it is rather small. Even during the whole period of its existence, T « 4.55 x 10 9 years, the 
sun has generated, via the three-body reaction (|3|), 

f^J « 0.835 x 10 33 erg , 

\ / era 

which is less than the energy radiated every second. This energy, however, has been stored 
inside the sun unless it was somehow transferred to the surface. In this connection the energy 
balance in the sun should, perhaps, be re-examined with inclusion not only the reaction (|3|) 
but also other nonradiative fusion reactions such as 



e + p + d — ► e + 3 He 
e + p + 7 Be — > e + 8 B , 
e + p + 7 Li — ► e + 8 Be 
etc. 



In conclusion, we can say that the three-body reaction (|3|) contributes to 7 Be production 
in the sun only about 0.01%, it perhaps plays a role in the energy balance in the sun, and 
definitely it was important at the early stages of the primordial nucleosynthesis. 
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TABLES 



TABLE I. Temperature dependence of the non-radiative capture rate, (TZ eTa )e : in units 
cm 6 mole _2 sec _1 . 



e (io 6 °k) 


(mi 


e (io 6 °k) 


(mi 


1 


4.653 


X 


10~ 55 


100 


2.656 x IO" 10 


2 


3.920 


X 


-i r\ — dl 

10 43 


110 


6.177 x 10~ 10 


3 


2.363 


X 


10 61 


120 


1.302 


1 n— 9 
X 10 


4 


1.039 


X 


10 66 


130 


2.531 


1 n — 9 

x 10 


5 


4.034 


X 


1 SI 

10 61 


140 


4.606 


i n— 9 
X 10 


6 


3.797 


X 


-i r» — 99 

10 za 


150 


7.927 


1 n — 9 
X 10 


7 


1.427 


X 


-i r» — 97 

10 li 


160 


1.301 


x 10~ 8 


8 


2.839 


X 


-i c\ — 9fi 

10 2b 


180 


3.119 


x 10 


9 


3.549 


X 


10 z& 


200 


6.605 


x 10~ 8 


10 


3.122 


X 


-i r* — 9/L 

10 24 


250 


2.944 


1 n — 7 

x 10 ' 


11 


2.087 


X 


10 2d 


300 


9.123 


1 n — 7 

x 10 ' 


12 


1.121 


X 


-i rv — 99 

10 22 


350 


2.236 


x 10~ 6 


13 


5.034 


X 


-i — 99 

10 22 


400 


4.661 


x 10~ 6 


14 


1.950 


X 


i r, — 91 

10 21 


450 


8.636 


x 10~ 6 


15 


6.671 


X 


i rv— 91 

10 Zi 


500 


1.464 


x 10~ 5 


16 


2.053 


X 


i r, — 9fl 

10 2U 


600 


3.469 


x 10~ 5 


18 


1.499 


X 


10 -19 


700 


6.853 


x 10~ 5 


20 


8.289 


X 


10 -19 


800 


1.195 


x 10~ 4 


25 


2.538 


X 


io- 17 


900 


1.901 


x 10~ 4 


30 


3.432 


X 


10 -16 


1000 


2.826 


x IO" 4 


40 


1.511 


X 


1Q -14 


1250 


6.170 


x IO" 4 


50 


2.208 


X 


io- 13 


1500 


1.104 


x 10~ 3 


60 


1.697 


X 


io- 12 


1750 


1.738 


x 10~ 3 


70 


8.606 


X 


io- 12 


2000 


2.507 


x IO" 3 


80 


3.273 


X 


io- 11 


2500 


4.380 


x 10~ 3 


90 


1.009 


X 


10 -io 


3000 


6.583 


x 10~ 3 
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TABLE II. Non-radiative reaction rate, (TZ eTa )e, its ratio to the corresponding radiative reac- 
tion rate, and the rate of the energy generation due to the non-radiative fusion, as functions of the 
radial position in the sun. 



R/Rq 


(R-eTa)e cm- 3 sec _1 






AE erj 


;/ (cm 3 sec) 


0.00 


7148 


0.294 


x 


10~ 3 


0.182 


x 10~ 13 


0.01 


7110 


0.288 


x 


10~ 3 


0.181 


x 10~ 13 


0.02 


6387 


0.282 


x 


10~ 3 


0.162 


x 10" 13 


0.03 


5416 


0.272 


x 


10~ 3 


0.138 


x 10~ 13 


0.04 


4139 


0.260 


x 


10~ 3 


0.105 


x 10~ 13 


0.05 


3202 


0.246 


x 


10~ 3 


0.814 


x 10~ 14 


0.06 


2194 


0.230 


X 


10~ 3 


0.558 


x 10~ 14 


0.07 


1407 


0.213 


X 


10~ 3 


0.358 


x 10~ 14 


0.08 


810 


0.196 


X 


10~ 3 


0.206 


x 10~ 14 


0.09 


515 


0.180 


X 


10~ 3 


0.131 


x IO" 14 


0.10 


270 


0.163 


X 


10~ 3 


0.688 


x 10~ 15 


0.15 


10 


0.962 


X 


io~ 4 


0.258 


x 10~ 16 


0.20 


0.241 


0.517 


X 


10~ 4 


0.612 


x 10~ 18 


0.30 


0.334 x 10~ 4 


0.123 


X 


10~ 4 


0.850 


x 10~ 22 


0.40 


0.126 x 10~ 9 


0.274 


X 


10~ 5 


0.321 


x 10~ 27 


0.50 


0.194 x 10" 14 


0.643 


X 


10~ 6 


0.493 


x 10~ 32 



TABLE III. Temperature dependence of the ratio of the three- to two-body reaction rates for 



the primordial nucleosynthesis. 



9 (10 6 °K) 




300 


0.181 X 10~ 8 


350 


0.404 x 10~ 7 


400 


0.428 x 10~ 6 


450 


0.275 x 10~ 5 


500 


0.124 x 10~ 4 


600 


0.124 x 10~ 3 


700 


0.671 x 10~ 3 


800 


0.244 x 10" 2 


900 


0.684 x 10~ 2 


1000 


0.159 x IO" 1 


1250 


0.766 x IO" 1 


1500 


0.231 


1750 


0.530 


2000 


1.02 


2500 


2.70 


3000 


5.44 
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FIGURES 

FIG. 1. Jacobi vectors in the configuration and momentum space. 




FIG. 2. Low-lying energy levels (in MeV) of 7 Be nucleus and the 3 He- 4 He threshold energy. 



1.5876 ( 3 He+ 4 He) 


6.73 (5/2-) 


4.57 (7/2-) 


0.4291 (1/2") 





7 Be (3/2-) g.s. 



15 



FIG. 3. Effective 3 He- 4 He potential for three different states corresponding to relative orbital 
angular momentum t = 0, 1 and total angular momentum J = 1/2, 3/2. 
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